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This research program investigates boundary-layer receptivity in the leading edge region for 
bodies with blunt leading edges. Receptivity theory provides the link between the unsteady distur- 
bance environment in the freestream and the initial amplitudes of instability waves in the boundary 
layer. This is a critical problem which must be addressed in order to develop more accurate pre- 
diction methods for boundary-layer transition. 

Previous analyses of leading-edge receptivity have been restricted to the case of a zero-thickness 
flat plate in symmetric flow at low Mach number. Singular perturbation techniques were utilized, 
based on a large Reynolds number for the unsteady motion (e = M ^) 1/6 <C 1). In this project 
we are extending the theory to include the effect of a parabolic leading edge, and also the effect of 
aerodynamic loading. The influence of nose bluntness enters through the parameter S = u;r n /t T co , 
where r n is the nose radius of the parabola. Our theory assumes that S = 0(1). 

Results for a symmetric base flow with a parallel acoustic disturbance are discussed in Ref. 1. 
Only the salient points are reproduced here. The boundary layer structure is illustrated in Fig. 1 . 
Parabolic coordinates 77 are used throughout. The streamwise and transverse coordinates £ and 77 
have been nondimensionalized by the convective wavelength, and by cHUo/uj, respectively. 

Unlike the flat-plate case, the mean boundary layer equation must be solved numerically. This 
was accomplished using a Keller Box method. The solution must then be matched to the large 
asymptotic limit, 
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Here F(g) is the Blasius (flat-plate) solution, g = gF* — F and / is the solution of a third-order 
ODE. The numerical constant A \ is determined by the condition that vorticity decays exponentially 
away from the body, but B\ can only be determined by comparison to the numerical solution. The 
extraction of B\ is illustrated in Fig. 2 . 

Far downstream (£ 1 ), the unsteady disturbance in the boundary layer is of the form 
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f= 1 

The functions xp, ip, are determined by the local properties of the geometry and the freestream dis- 
turbance in the far-downstream region. For a n a c ous tic freestream disturbance, rp corresponds to a 
Stokes shear w r ave. The asymptotic eigensolutions depend only on the downstream geometry and 
are independent of theTreestream, disturbance. However, the coefficients C,- are determined by the 
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receptivity process which occurs in the region £ — 0(1). Specifically, the first asymptotic eigenso- 
lution u > i matches onto the Tollmien-Schlichting wave which becomes unstable farther downstream 
in the boundary layer. Thus it is through Ci(S) that the freestream disturbance determines the 
magnitude of the unstable disturbance, and hence C x is called the receptivity coefficient. 

The asymptotic eigensolutions 0,- decay exponentially with downstream distance and they are 
also inverse ordered, i.e., the ratio ipi/ip i+ x becomes exponentially small for £ > 1. Thus, it is 

difficult to obtain C x from a direct numerical solution of the Linearised Unsteady Boundary Layer 
Equation (LUBLE). This problem is circumvented by extending the analysis and the numerical 
solution to the complex £ plane and choosing arg(f) such that p x grows exponentially. In this 
region of the complex plane the asymptotic eigensolutions assume a normal ordering, with rp x 
dominating both $ and the higher terms in the above series. The Receptivity Coefficient C 1 (5) can 
then be obtained by comparing the numerical results with the asymptotic form of tp x . Essentially, 
the known analytical form of the asymptotic eigensolution is utilized to extrapolate the numerical 
results to £ — oo. For £ > 1, the first asymptotic eigensolution has the form 

0i ~ (i + b,r n 774 + hr 1 log 2 <e + hr 1 ^ + & 4 c‘ + . . o cK*?)e rf ° (3) 
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The fact that there are several terms in (3) which are similar in magnitude makes the extrapolation 
difficult, and numerical integration must be continued up to very large values of £ to obtain accurate 
estimates of the Receptivity Coefficient C\. 

The variation of C x with S for an acoustic disturbance parallel to the body axis was given in 
Ref. 1. A dramatic decrease in receptivity with increasing nose radius was reported. At S = 0.3 the 
Receptivity Coefficient C\ drops to 10% of that for a flat plate. The conclusion to be drawn is that 
a favorable pressure gradient not only shifts the point of neutral stability further downstream, but 
also inhibits the receptivity process. An asymptotic theory for small S has now been completed. In 
this case the region close to the nose is quasi-steady and in the region £ = 0(1), which dominates 
the receptivity process, the body surface approaches a flat plate. This analysis predicts an increase 
in receptivity for small S , 

\Ci(S)\~ |<7 1 (0)|(1 + 3.135). (4) 


Comparison with numerical results (Fig. 3) shows that this is indeed true, but also reveals that 
the 5 < 1 theory is useful only for very small $ and that the maximum increase in \C X \ is only of 
the order of 3%. The reason behind this increase is not fully understood. However, for all practical 
purposes this very local increase can be ignored. 

The Receptivity Coefficient has also been calculated for different types of freestream distur- 
bances. In carrying out these calculations we take advantage of the fact that the evolution of the 
unsteady disturbance (2) with downstream distance in the complex £ plane occurs in two distinct 
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stages. During the first stage the unsteady disturbance becomes dominated by the single term 
Ci (5)^i (>?.£, S), which grows exponentially relative to the other terms in (2). Due to the exponen- 
tial dependence only moderately large values of £ are necessary for this first stage to be completed. 
The second stage corresponds to the evolution of the function t/q according to the asymptotic form 
(3). The higher-order terms in (3a) are algebraic and therefore become unimportant only for much 
larger values of £. In extracting the Receptivity Coefficient for an acoustic disturbance parallel 
to the body axis, it was necessary to continue the numerical computation out to distances where 
(3) is dominated by the first few terms. However, since this computation provides information 
regarding the evolution of 0 l9 in the computations for other types of freestream disturbances it is 
only necessary to extend the computations to the point where the first stage has been completed. 

For an incoming acoustic disturbance at angle 9 to the axis of symmetry, the Receptivity 
Coefficient is given by 

Ci (5) = cos 9C a (5) + a 1/2 sin 6C b (S) (5) 

where a is the airfoil half-chord, non-dimensionalized by the convective wavelength. This is valid for 
a > e“ 2 , ensuring that transition occurs towards the leading edge. The magnitudes and arguments 
of C a and C b are plotted in Figs. 4 and 5. It is seen that the receptivity due to the asymmetric 
component of the disturbance is much larger than that due to the symmetric component, but it also 
reduces faster with S than the symmetric component. This is illustrated more clearly in Figures 6 
and 7, where the total receptivity is illustrated for a realistic non-dimensional chord length, a = 225. 
In Fig. 6 the variation [Ci| with incidence angle is illustrated for various values of nose radius. In 
Fig. 7 the modulus of the receptivity is plotted against S for various acoustic incidence angles. For 
small values of 5, receptivity levels rise rapidly with increasing angle of incidence of the acoustic 
wave. However at larger values of 5, this increase is less marked. These results are for an airfoil 
with no induced circulation due to the oblique disturbance. For a body with sharp trailing edge, 
an additional contribution to the slip velocity due to enforcement of the unsteady Kutta condition 
must be included. This will be calculated in the near future. 

Work is now proceeding to determine the receptivity due to a convected vorticity gust. In 
contrast to the flat-plate case, the vorticity is distorted by gradients in the inviscid flow field near 
the leading edge. The slip velocity on the body due to the gust is obtained using Rapid Distortion 
Theory. 

* 

We are also extending our theory to incorporate the influence of mean aerodynamic loading on 
receptivity. This effect has been ignored in previous studies because loading is incompatible with 
the zero-thickness, flat-plate model. Aerodynamic loading produces an adverse pressure gradient on 
the upper surface, and we anticipate that this may lead to a substantial increase in the Receptivity 
Coefficient. 

The asymptotic analysis of the steady and unsteady boundary layer flow in the nose region, 
for an airfoil at non-zero angle of attack has now been completed. This will be published in detail 
elsewhere. Here we give a few of the more important results. For an angle of attack a, the steady 
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boundary layer solution well downstream of the nose is given by 
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where // = sin a(a/S) 1 ^ 2 and w = cosa£ + //, and f,g,p,q,r are functions of rj given by a set of 
third-order ODEs. The stagnation point is given by w = 0. 

As in the symmetric case, A x is determined by the condition of exponential decay of vorticity, 
but now depends on p. B x {p) is determined by comparing numerical solutions to the asymptotic 
solution. This has been carried out and B x obtained for various values of fi. All the terms included in 
the expansion of ^ are necessary in order to calculate the eigenfunction for the unsteady disturbance 
in the boundary layer. However, the need to include so many terms in the expansion naturally 
makes the analysis considerably more involved. Finally the asymptotic eigenfunction is found to 
be 
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and a,6,c,d are known functions of /i. The computer code used for the symmetric case is currently 
being modified to solve this modified LUBLE. The first results are expected to be available in the 
next month. For acoustic disturbances in the freestream, the coefficient of receptivity will now be 
a function of three parameters, angle of attack, nose radius and acoustic angle of incidence. 

In summary, our research to date has determined the influence of nose radius on leading-edge 
receptivity for the case of symmetric flow past a parabolic leading edge. Results have been obtained 
for a variety of acoustic freestream disturbances, and calculations for the case of vortical freestream 
disturbances are presently being carried out. In addition, the asymptotic theory has been extended 
to include the influence of mean aerodynamic loading in the leading-edge region, and the related 
computer programs are presently being developed. This work will provide answers to fundamental 
questions concerning leading-edge receptivity, i.e,, the influence of the airfoil nose radius and mean 
aerodynamic loading. These issues must be understood in order to apply receptivity theory to 
realistic geometries. 
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Figure Captions 


Figure L A schematic illustration of the boundary layer structure 
for a parabolic body in a symmetric mean flow. The three decks 
in the Orr-Sommerfeld region are (i) the viscous wall layer; (ii) the 
main inviscid layer; and (iii) the outer irrotational layer. 

Figure 2. Plot of G(Z) = (log Zy l Z(\4>" - F"]}F"\= 0 against 
(logZ) -1 , where Z = w 2 . Comparison with the large w asymptotic 
form yields & 2.08. 

Figure 3. Comparison of asymptotic and numerical results for very 
small S . Numerical results are marked with squares. 

Figure 4. Variation of the amplitude of the Receptivity Coefficients 
as a function of Strouhal number, S. C a and C& are the Receptiv- 
ity Coefficients for the symmetric and antisymmetric components of 
freestream disturbance respectively. 

Figure 5. Variation of the phase of the Receptivity Coefficients as a 
function of Strouhal number, $. 

Figure 6. Variation of |C| with acoustic incidence angle 9 , for non- 
dimensional airfoil chord length a = 225. 

Figure 7. Variation of \C\ with S for various acoustic incidence 
angles, with the same airfoil chord length as Figure 6. 
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Figure 1. A schematic illustration of the boundary layer structure 
for a parabolic body in a symmetric mean flow. The three decks 
m the Orr-Sommerfeld region are (i) the viscous wall layer; (ii) the 
main inviscid layer; and (iii) the outer irrotational layer. 





Figure 4. Variation of the amplitude of the Receptivity Coefficients 
as a function of Strouhal number, 5. C a and C are the Receptiv- 
ity Coefficients for the symmetric and antisymmetric components of 
freestream disturbance respectively. 



Figure 5. Variation of the phase of the Receptivity Coefficients as a 
function of Strouhal number, S. 



Figure 6. Variation of |C'| with acoustic incidence angle 6. for non- 
dimensional airfoil chord length a = 225. 
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Figure 7. Variation of |C| with 5 for various acoustic incidence 
angles, with the same airfoil chord length as Figure 6. 
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